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Kicked Harper operators and on-resonance double kicked rotor operators model quantum systems 
whose semiclassical limits exhibit chaotic dynamics. Recent computational studies indicate a striking 
resemblance between the spectrums of these operators. In this paper we apply C*-algebra methods 
to explain this resemblance. We show that each pair of corresponding operators belong to a common 
rotation C* -algebra B a , prove that their spectrums are equal if a is irrational, and prove that the 
Hausdorff distance between their spectrums converges to zero as q increases if a = p/q with p and q 
coprime integers. Moreover, we show that corresponding operators in B a are homomorphic images 
of mother operators in the universal rotation C*-algebra A a that are unitarily equivalent and hence 
have identical spectrums. These results extend analogous results for almost Mathieu operators. We 
also utilize the C*-algebraic framework to develop efficient algorithms to compute the spectrums 
of these mother operators for rational a and present preliminary numerical results that support the 
conjecture that their spectrums are Cantor sets if a is irrational. This conjecture for almost Mathieu 
operators, called the Ten Martini Problem, was recently proved after intensive efforts over several 
decades. This proof for the almost Mathieu operators utilized transfer matrix methods, which do 
not exist for the kicked operators. We outline a strategy, based on a special property of loop groups 
of semisimple Lie groups, to prove this conjecture for the kicked operators. 
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I. INTRODUCTION 

The families of kicked Harper and on-resonance double kicked rotor operators that we discuss in this paper 
arose comparatively recently in the field of quantum chaos and most of the knowledge about them has been 
acquired through numerical computation. They are related to another family of operators that arose earlier, that 
is simpler, and whose properties are far better understood theoretically. We start by reviewing this family of operators. 

We let C, R, Z, and N = {1,2,3,4, ...} denote the complex, real, integer, and natural numbers. Furthermore, we let 
T = R/ Z denote the circle group parameterized by the interval [0, 1), and T c denote the circle group parameterized 
by the numbers on the complex unit circle. The family of almost Mathieu self-adjoint operators are represented on 
L 2 (T). The representation with respect to the standard orthonormal basis {£, n (x) = exp(i2irnx) : n S Z } is 

H(a, \,6) £ n = £ n+1 +£ n - 1 + 2\ cos [2>ir(na + 6)] £„, n £ Z. (1) 

Here the coupling constant A e R\{0}, the frequency a e [0, 1), and the phase 9 G [0, 1). These (also called Harper) 



operators were introduced in 1955 by Harper [27j to explain the effect of magnetic fields on conduction bands of a 
metal. The spectrums of these operators have been studied for several decades and are well understood. In 1964 
Azbel conjectured that if a is irrational then the spectrum of H(a,X,6) is a Cantor set. In 1976 Hofstadter [3], 
observing that the spectrum of H (a, A, 9) is independent of 9 if and only if a is irrational, argued that "there can be 
no physical effect stemming from the irrationality of some parameter" and proposed to study instead the set 

S(a,X)= [j spectrum [ff (a, A, 0)]. (2) 
ee[o,i) 

He computed S(a, 1) numerically for a fine grid of rational values of a € (0, 1) and studied the graph of this set- valued 
function of a (a subset of the rectangle [—4, 4] x (0, 1)). He asserted that the graph had a butterfly pattern, discussed 
the relationship between the recursive clustering structure of the graph and the continued fraction expansion of a. In 
1981 the problem of proving Azbel's conjecture was named the Ten Martini Problem by Barry Simon after an offer 
by Mark Kac [47]. In 1982 Simon and Bellisard 0] showed that if a = p/q where p and q are coprime, then S(a, A) 
consists of q disjoint closed intervals if q is odd and q — 1 disjoint closed intervals if q is even, and they used this result 
to prove that for a Gs set of pairs (a, A) in R 2 , S(a, A) is a Cantor set. In 1993 Choi, Elliot and Yui 16| used the 



noncommutative binomial theorem to compute sharper estimates for the spectral gaps when a is rational and used 



these estimates to prove that S(a, A) is a Cantor set whenever a is a Liouville number. In 1993 Last 32| proved that 
for almost all a 6 (0, 1) the Lebesque measure of S(a, A) equals 4 1 1 — |A| | for all A € R, and hence that S(a, 1) is a 
Cantor set. In 1987 Sinai [3| used KAM theory to prove that S(a, A) is a Cantor set for almost all a and sufficiently 
small (or large) |A|. The Ten Martini Problem was recently proved by Avila and Jitomirskaya Closely related to 
the Harper operators, we define the family of unitary Harper operators 

Uh{k, a, A, 9) — exp [ — i k H(a, A, 9) ] (3) 

where n £ R and a, A, 9 are as in Equation [TJ We will show in Corollary [5] that if a is irrational, then the spectrum 
of Uh is independent of 9 and is a Cantor set. 

In 1979 Casati, Chirikov, Ford and Izrailev [3| and Berry, Balazs, Taber and Voros [8] initiated the study of 
families of unitary operators that describe the time evolution of quantum systems whose classical limits exhibit 
chaotic dynamics. Developments in the field of quantum chaos [l5j led to the so-called kicked Harper operators that 
admit the following representations on L 2 (T) : 

Uk{i(K, a, A, 9) — exp — i2n cos(27rx) exp — i2n A cos ^ — ia— — h 2tt9^ . (4) 
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We observe that t/;yj(K, a, A, 6) — > Uh(k, a, A, 0) as |«| — > 0. In 1990 Leboeuf, Kurchan, Feingold and Arovas [34 1 
studied a quantum version of a classical kicked Harper map whose phase space is a (compact) torus and that exhibits 
chaotic dynamics for sufficiently large values of certain parameters. Their quantum kicked Harper operators arise 
through canonical quantization of this classical Harper map and have the form Ukni^, ot, A), (see Equation 1281 below) 
and is represented on the Hilbert space L 2 (T 2 ). To quantize the system on the torus they set a to be the area of classical 
phase space in units of Planck's constant and they choose a = 1/q with q an integer. However, their quantum operators 
admit the family (parameterized by 8) of representations on L 2 (T) that have the form Ukn(K, 1/q, A, 8). In 1991 Geisel, 
Ketzmerick and Petschel [22J studied the dependence on k of the spectra of the operators Ukii( K , a, 1), where a ranges 
over the interval [0, 1). These spectra are the same as the union over 8 of the spectrums of a, A, 8). They also 

remark that 9 can be thought of as being proportional to a quasimomentum which arises in analogy to Bloch's theorem. 
They observe that for small |«| the collection of spectrums for rational a G [0, 1) has the Hofstadter butterfly pattern 
for small |k|, and that as \n\ increases the spectral bands close up. They also compute the effect of K on other 
quantities such as diffusion and localization. In 1992 Artuso, Casati, and Shepelyansky Q studied a quantum version 
of another classical kicked Harper map whose phase space is a (noncompact) cylinder. Their operators have the form 
UkH^K, a, A, 0). This quantization evidently differs from the canonical quantization of Leboeuf et at insofar as the 
state space of the quantum system equals L 2 (T), the quasimomentum parameter is fixed at 9 — 0, and there are no 
restrictions on a. The quantum version with 9 = was later widely used, though a more general quantization gives 



the operator UkH («, cc, A, 9) [12| . The fact that a classical kicked Harper map can be quantized in many ways was 



probably best summarized by Guarneri and Borgonovi 25 1 



Another family of operators we study here is the so-called on-resonance double kicked rotor operators. They admit 
the following representations on L 2 (T): 

Uordkr(n, a, A, 9) — exp — %2k cos(27tx) T(a) exp — i2n A cos (2n(x + 8)) T(a)~ , (5) 

where T(a) = exp ~ix^) . The operators U or dkr{^, ct, A, 0) were introduced by Gong and Wang in 24|, where 
they compared their properties with those of operators a, A, 0). In [52| they discussed how to implement the 



operators U„ r dkr(K, ol, A, 8) for 8^0 and related 8 to a symmetry breaking condition in a Hamiltonian ratchet model. 



In 51| and 5fJ they and Mouritzen discussed how the operators U or dkr{^, a*, A, 0) can be realized experimentally, and 
demonstrated numerically a striking resemblance between the spectrums of i7fcg(l/2, a, 1, 0) and [/ or dfer(l/2, a, 1, 0) 
for a large set of rational a = p/q with p and q coprime. In particular, Figure 2 in [5l| shows that the spectrums for 
a = 13/41 are different, but very nearly equal. 

One main objective of this work is to mathematically explain this resemblance between on-resonance double kicked 
rotor operators and the kicked Harper operators, by proving that for all values of (k, X, 81,82) the spectrums of 
U or dkr(n, a, A, 8\) and J7feH(/t, a, A, 82) are equal whenever a is irrational and that their Hausdorff distance approaches 
zero as q increases for rational a = p/q. Because Wang and Gong [5l| also showed that the dynamics of the on-resonance 
double kicked rotor model dramatically differs from that of the kicked Harper model for irrational a, our proof of 
the spectral equivalence, in addition to mathematical interest, shall motivate new theoretical problems, such as the 
characterization of (generalized) eigenfunctions of the operators and of density of states or of the degeneracy of the 
spectrums, which ultimately accounts for the dynamical differences observed by Wang and Gong. 

Our proof is based on C* -algebra methods. The proof itself also leads to efficient algorithms for numerical stud- 
ies of the above-mentioned unitary operators. Furthermore, we discuss the observation that the spectral graph of 
U or dkr{n, a, A) has a butterfly pattern similar to that previously observed for Ukui^i a, A). Though rigorous spectral 
characterization of UkH (k, a, A) and U or dkr{i^, ol, A) remains an open problem, we outline a strategy to prove the con- 
jecture that the spectrum of their mother operators Ukn{n, ci, A), or identically U or dkr{^, ot, A), is indeed a Cantor set 
for irrational a. These mother operators, denoted with a tilde, live in the universal rotation C*-algebra A a and the 
corresponding operators Ukii(K, A) and U or dkr(K, a, A), living in the rotation algebra B a , are homomorphic images 
of them. 



4 



Our paper is arranged as follows: In section pi] we derive results concerning the spectrums of H, Uh, UkH and 
Uordkr, in particular the main result that for irrational a, the spectrums of U or dkr(K, a, X,8) and f7fcjj(/c, a, X,9) are 
identical. In section HTT1 we present efficient algorithms for numerically calculating the spectra of U or dkr(K, cc, A, 9) 
and UkH {k, a, A, 9) for rational a — p/q. The treatment in both of these sections HI1 and IlIII rely on results given 
in appendices [B] and O In section IIV( we discusses ideas for future research. Appendix [X] briefly reviews the 
considerations and difficulties involved in physically realizing the operators above. Appendix [5] summarizes results 
in spectral theory and the theory of C*-algebras used in the paper. Appendix [Cl further builds on these results and 
summarizes results regarding rotation C* -algebras used in the paper. 



II. DERIVATIONS 



For r > we define the disc D(r) = { z G C : \z\ < r }. For S C C we let H(S) denote the metric space whose 
points are compact subsets of S and that is equipped with the Hausdorff metric 



d(X, Y) = max { maxmin \x — y\ , maxmin \x — y\ | , X,Y £ Tl(S). 

1 x£X y£Y y£Y x£X 



(6) 



This metric space is called the Hyperspace of S. A standard result in topology ensures that the hyperspace H(S) 
is compact whenever S is compact ( 19|, p. 205, p. 253), ( 36J, p. 279). For every Hilbert space H we denote the set 
of bounded operators on H by B(H), its subset of operators that are normal by B n (H), unitary by £> U (H), and 
self-adjoint by B a (H). We observe that the real part of A 6 B(H), given by 5i(A) = ^(A + A*), is self-adjoint. If 
A e S a (H) then exp(iyl) g £> U (H) and Lemma 3.3 in Q implies that 



||exp(a)-exp(iB)|| < \\A-B\\, A,B e B a (H). 

We define a : B(H) -> H(C) by 

cr(B) = spectrum(£?) = {/i£C : /i/ - B does not have an inverse in £>(H) }, 



B e 6(H). 



(7) 



(8) 



As noted in the Examples following Definition QTJ B(H) is a C*-algebra, and Lemma [T71 as noted in the comment fol- 
lowing it, implies that for every B e B(H), cr{B) = crg( H )(i3). Therefore, Lemma [T^] implies that a(B) £ H\T)(\ \B\ |)) 
for every B G 6(H). Furthermore, a(B a (H)) C H(R) and cr(B„(H)) C H(T C ). If dim(H) < oo and A S B(H) then A 
can be represented by a matrix and <r(A) is the set of eigenvalues of this matrix. This is not the case when H is infinite 
dimensional. For example, if H = L 2 ([0, 1]) and (Af)(x) — xf(x) then a(A) — [0, 1] and A has no eigenvalues since 
it has no eigenvectors. Although a is not continuous ([3l|, Example 3.8), Proposition [T91 shows that the restriction of 
a to B n (H) is continuous since 



d[a(A),a(B)} < \\A-B\\, A,Be B„(H). 



(9) 



We observe that since the operators H, UkH, and U or dkr are normal, Equations and © imply that their spectrums 
depend continuously on the parameters k, A, 9. In particular 



d 



a(H(a,X,6x)), a(H(a,X,6 2 )) < 2] A] | sin[7r(^ - 6- 



a(U H (K,a,X,e 1 )),a(U H (K,a,X,e 2 )) < 2\k X\ \ sin^x - 9 2 )] 
a(U kH (K,a,X,e 1 )),a(U k H(K,a,X,e 2 ))] < 2\k X\ \ sin[7r(^ - 9 2 )] 
cr(/7 or(i fcr(K,a, A, 6<i)), a(U or dkr(K,a,X,9 2 )) < 2\k X\ \ sin[7r(0i - 6 2 )] 



(10) 

(11) 
(12) 
(13) 
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Remark 1 We observe that if ot\ ^ Z, ct2 Z, a\ — cti ^ Z, one? ai + o?2 ^ Z. i/ien 



|| H(ai, A, 0) — H(a2, A, 0) || = sup 2A sin [7rn(a!i + 012) + 2n9] sin [7rn(ai — 0:2)] 

Trais shows that the dependence on a is not necessarily continuous. 
We define M : L°°(T) -> 6(i 2 (T)), i? : T -> tf(L 2 (T)), by 

(M( 5 ) /)(*) - ff (x)/(a:), 9 £ L°°(T), / e i 2 (T), ieT, 
(i?(a) /) (x) - f(x + a), a £ T, / e i 2 (T), .x e T, 



> #|A| 



(14) 



(15) 
(16) 



and let B a C £>(L 2 (T)) denote the rotation C* -algebra generated by the operators M(£i) and R(a) and their adjoints. 
Then (R(a), M(fi)) is a frame with parameter a for i3 Q , see Definition (53] 



Proposition 2 77ie operators H(a, \,9), Uh(k,c<,\,0), £/&#(«;, a, A, and U or dkr(K, ct. A, 0) belong to B a since 

2R[M(£ 1 )]+2AR[fc(0)ii(a)] l (17) 

exp [ - i2«5Jt(M(fi)) - i2 K A5R(Ci ((9) fl(a))] , (18) 

exp [ - i2u3fi(M(fi))] exp [ - i2nm(^ {9) R{a))) , (19) 

exp [-i2«3fi(M(fi))] exp [ - «2kAK(^i(0 + a/2) M(fi) (20) 



H(a,\,6) 
U H {K,a,\,9) 
U kH {K,a,\,9) 
U or dkr(n, a, A, 0) 



Proof. The first equation follows from 23? [M(£i)] £„ = £„+i + fr»-i and 2A3? [6(6*) #(a)] = 2Acos [27r(na + 9)] £„. 
The second equation follows from (— ia-^ + 2n9) £ n = 2n(na + 6)£ n . We observe that T(a) £ n = exp {iiran 2 ) £„ and 
hence the third equation follows from T(a) M(£i) T(a) _1 ^„ = exp(— i7ra) i?(a) M(£i) £ n . ■ 

We say that A, B e ^(H) are unitarily equivalent (and write A = B) if there exists an O £ B U (H) such 
that A = OBO^ 1 . Clearly A = B implies that a (A) — cr(B). We say that A and B are approximately unitarily 
equivalent (and write A = a B) if there exists a sequence Ofc € B U (H) such that limfc^oo ||A — OkBO^ 1 ]] = 0. 
Although approximate unitary equivalence is weaker than unitary equivalence [18|, if A, B £ B n (YL) and A = a B 
then <j(A) = <j{B). This follows since Proposition implies that a : B n (H) — > H(C) is continuous. 



Proposition 3 If a —p/q is rational and p and q are coprime, then for all #i,#2 € [0,1), 

a(H(a,\,e 1 )),a(H{a,X,e 2 ))] < 2|A| | sin[7r/(2g)]|, 
o-{Uh{k,ol,\,6i)), o-(U H {K,a,\e 2 ))\ < 2\k A| | sin[7r/(2g)] 
CT(cr fcH («,a s A,fli)),<r(l7 fcH (K,Qt ) A,ea))] < 2|« A| | sin[7r/(2g)] 
<r(Uordkr(K,oi, X,6i)), cr{U or dkr{K,a, A,0 2 )) < 2|kA| I sin[7r/(2g)] 



(21) 
(22) 
(23) 
(24) 



//a is irrational then the spectrums of these operators are independent of 8 and therefore the distances above are zero. 

Proof. The unitary equivalence M(£„) R(a) M (Cn) _1 = £i(~ na ) R( a ) implies that H(a,X,9) = H(a, A, 9 — no), 
Uii(K,a,X,9) = Uh(k,oi, A, 9 — na), £/&#-(«, a, A, 9) = Ukui^, a, A, — not) and that U or dkr{^iCt,\,9) = 
U rdkr(K, a, A, — na). Therefore, Inequality 1101 implies Inequality 1211 Inequality 1111 implies Inequality 1221 Inequal- 
ity [12] implies Inequality [23l and Inequality [13] implies Inequality [24] The last statement follows since if a is irrational 
then aZ c T is dense and hence the unitary equivalences above become almost unitary equivalences. ■ 

Corollary 4 If a is irrational, then the spectrum of Uh{k, a, A, 9) is independent of 9 and is a Cantor set. 

Proof. Proposition [3] implies the first assertion. Then Lemma [TBI implies that cr(J7^r(K, a, A, 0)) is a Cantor set since 
it equals the image of the Cantor set o~(H(k, a, A, 9)) under the map t — > exp(— int). M 
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Following Proposition [211 for a G [0,1) we let (U a ,V a ) be a frame with parameter a for the universal rota- 
tion C*-algebra A a and we construct operators, called mother operators, in A a by 



H(a,X) 
U H (n,a,X) 

U kH (K,a,X) 

Uordkr(K, a, A) 



exp(— iira) U a V a , 
23t(V a ) +2Xdi(U a ), 



exp 
exp 
exp 



-i2K$t(V a ) 
-i2n%t(V a ) 



exp 
exp 



-i2nXtl(U a ) 
-i2nXVi(W a ) 



(25) 
(26) 
(27) 

(28) 

(29) 



Theorem 5 The families of operators defined by the equations above satisfy 



% e (H(a,\)) = H(a,X,6), 
irg(U H (K,a,\)) = U H {n,a,X,9), 
ir e (U kH (K,a,X)) = U kH {n,a,\,9), 

^e(U or dkr(K,a,X)) = U or dkr(K,a,X,9), 

where -kq 1 9 € T € [0, 1), is the family of homomorphisms constructed in Lemma 

a(H(a,X)) = 



(30) 

(31) 
(32) 
(33) 



and their spectrums satisfy 



o-(U H (K,a,\)) 



|J a(H(a,X,0)), 
ee[o,|) 

|J a(U H (K,a,\,6)), 
8e[o,i) 

|J o-(U kH (n,a,\,6)), 
ee[o,i) 

[J 0-{U O rdkr(K, a, A, 9)) 

ee[o,i) 



and are continuous functions of n, a, and A since 

d ' 

d 

d ' 



a(H(a!,\)), a(H(a 2 ,X)) < 36 V67r|A(a 2 - ai )| , 

a(U H (k, ai, A)), a(U H (n, a 2 , A)) < 36 \Z67t|kA(q?2 - ai)| 

o-(U k H(K, at, A)), a(U kH (K, a 2 , A)) < 36 v / 67t|kA(c^ - Oij] 

o-(U ordkr {K, ax, A)), o-(U ordkr (K, a 2 , A)) < 36 ^67r|KA(of 2 - ai)| 

Furthermore, there exists an L E B{L 2 {T 2 )j such that 

U kH (K, a, A) = LU kH (n,a,X) L^ 1 = U ordkr {n>,OL,X) 
and hence a(U k H{n,a, X)) = a(U or d kr (K,ct, X)) . 



(34) 
(35) 
(36) 
(37) 

(38) 
(39) 
(40) 
(41) 

(42) 



Proof. Proposition [2] and Lemma [30] implies that the homomorphism irg maps the operators defined by Equations 
[26H291 onto the operators on the right side of Equations I301I331 The proof of Proposition [3] shows that each of the 
operators that appear on the right in Equations I34hl37l is unitarily equivalent to itself with 9 replaced by 9 + 1/q. 
Therefore, Proposition |2"01 implies Equations |3"4"H3"T1 Let /i > and let Wj : A aj — > B(H) for j = 1,2 be the injective 
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homomorphisms described in Proposition [32l Define operators Hj € B(H), j = 1,2 by Hj = nj (H(ctj, A)) . Since iXj is 
injective a[H(aj, A)) = cr(Hj) and hence Inequality O implies that d[cr(i/(ai, A)), cr(-ff(a:2, A))] = d[a(H 2 ), a{H\)\ < 
\\H 2 — Hi\\. The triangle inequality together with Inequalities IC15I implies that 

~M ~ 54 

ff 2 -ffi < 2 3?tt 2 (V aa )-M7ri F ai +2 A 3?tt 2 E/„ 2 -^(t/^) < — + 18 A 

M 

Choosing = •v/3/|A] minimizes the right side of the inequality above and gives Inequality Inequalities |3"5H4"T1 are 
obtained using the same procedure. Proposition |2"51 implies Equation 1421 and completes the proof. ■ 



III. ALGORITHMS AND NUMERICAL EXAMPLES 

Assume that a — p/ q, where p € Z, q G N, with p and q coprime. We will derive algorithms that accurately estimate 
the spectrums of the unitary operators Uh(h, a, A, 6*), £/&#(«;, a, A, 6*) and U or dkr(K,a, A, 6*) and of their mother 
unitary operators Uh(k, ct, A), t//c#(/-t, a, A) and U or dkr{i^, ct, A). The algorithms work by computing eigenvalues 
of matrices in A4 q , the set of q x q matrices. The matrices are parameterized by (x,9) E [0, 1/q) 2 and the es- 
timates are obtained by forming the union of the sets of eigenvalues over an equal-spaced NxN grid of values of (x, 9). 

Let (j = exp(i27r/q) and let F, C, and D be the q x q matrices defined by Equations IC21 For x € [0,1) con- 
struct the homomorphism p{x) S Hom(So,,A^ g ) so that p(x)(M((,i)^ = £i(x)D and p(x)(R a ) — C p . Proposition [27] 
ensures the existence of this homomorphism. For any integer k and y € [0, 1) define the function c(y) = cos(2iry) and 
define the q x q matrix 



G(k,y) 



c{y) 



( , k(q-l)^ 



The following result provides the basis for the algorithms for the unitary Harper and unitary kicked Harper operators. 



Theorem 6 For a = p/q, [0, 1) and K,AeRuc define the matrices 

M UH (K,a,X,x,9) = exp \-i2n{G{l,x) + XFG(p, 9) F^ 1 ) 
and 

M UkH (k, a, A, x, 9) = exp [ - i2n G(l, x)] F exp [ — i2n\ G(p, 9)] F" 1 

-i2nc(x) ] |~ -i2 K \c(B) 



F- 1 



-i2n c[ x+- 



F 



(43) 



(44) 



-i2 K \c(e+ p{q - 1) ) 



' Vl 


v 2 ■ 


■ V q -1 


Vq 


Vq 


Vl ■ 


■ 


Vq-1 


V3 


v 4 ■ 


• Vi 


V2 


. V2 


v 3 ■ 


■ V q 


Vl . 



-i2n\c{9) 



-i2fcAc(e+£) 



-i2 K Ac(e+-^^£) 



-i2k\c 9+i 
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where 



[vi v 2 v 3 



■ F 



-i2«c(x) -i2Kc(a;+i) -i2 K c(x+|) -i2ftA c(a;+i2_ii) 



If a = p/q, 9 £ [0, 1) and k, A G R, i/ien 



a(U H (K,a,X,9)) = (J <t(M Uh (k, a, A, x, 6*)) , 



(45) 



a(U kH (K,a,X,6)) = (J cr(Mu hH (k, a, A, 0)) . 

//a uniformly spaced grid of N 6 N values of x is used to compute an estimate o~N\U k H (k, a, A, 0)) t/ien 

27t|k| 



(46) 



cr N (U H {K,a,\,9)), <j(U h (k, a, A,6»)) 



< 



(47) 



a N (U kH (K, a, A, 0)) , a(U kH {n, a, A, 0)) 



< 



2tt|k| 



//a = p/q and k, A € R, i/ien 



and 



o"(?7ff(«»a,A)) = [J [J («;,«, A, £,0)), 

xe[o,i) 0e[o,i) 



o-(U k H(K,az, A)) = (J [J a(M UkH (k, a, A, x, 0)) . 

*e[o A) ee[o,i) 



(48) 



(49) 



(50) 



Furthermore, if a uniformly spaced two-dimensional grid of N x N values of (x,9) £ [0, i) 2 is used £o compute an 
estimate aNxN(UkH{n,ct: X)) then 



and 



o-nxn(U h (k,o>,\)), a(U H (n,a,X)) 



o- NxN (U k H(K,a,X)), a(U kH (K,a,X)) 



< 



< 



2ttH(1 + |A|) 
Nq 

27r|«|(l + |A|) 



Nq 



(51) 



(52) 



Proof. Since p(x)(R a ) = C p = F D p F^ 1 it follows that p(x)(U H (K,a, X,6)) = M UH (n,a,X,9) and 
p(x)(U k H(K, a, X, 0)) = M UkH (n,a, X,9). Therefore, since C G(l,x)C~ 1 = G(l,x + l/q), Proposition [2171 implies 
Equations US] and 05J Furthermore, Equations [35J and [35] imply Equations 25] and [55] Inequalities 23 and 2E1 are 
derived from Equation [Qj using Equation [Jj to show that 

\\M(K,a, X,xi,9) - M(n,a, A,x 2 ,0)|| < 2\n\\\G(l,x 1 ) - G(l,x 2 )\\ < 4tt\k(x 2 - xx)\, 

and the fact that the distance between every x £ [0, 1/q) and a uniform grid with N points in [0, 1/q) is < l/(2Nq). 
Inequalities 1511 and 1521 are derived likewise from Equation [9] using Equation [Jj to show that 

||M(K,a,A,a:i,0i)-M(K,a,A,a:a,02)|| < 2|/c| ||G(1, an) - G(l,a; 2 )| | + 2|kA| | \G{p, Ox) - G(p, 2 )| | 

< 2tt\k(x 1 -x 2 )\+2tt\kX(9 1 -9 2 )\. 
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This concludes the proof. ■ 



We now consider on-resonance double kicked rotor operators. Clearly, 



p(x)(U rdkr(K, a, A, 0)) 



exp 



i2kG(1,x) 



exp 



*2kA ft( + a/2) a (x) D C p ) 



(53) 



In order to develop an efficient algorithm we need to diagonalize the matrix DC P . Let v € C be an eigenvalue of 
DC P . Then there exists a (nonzero) eigenvector u — \u\ U2 
of u modulo q, compute 



u q ] T such that DC p u = vu. We consider the indices 



u p+ i 



Uqp+l 



1/1*1, 



UJ H V Ui, 



"3p+l = W P 2 ' P V 3 Ui, 
u ip+1 = u;- p - 2p - 3p v 4 Ul , 



UJ 



-p—2p—3p- 



Ui , 



(54) 



and then observe that DC p u = vu if and only if uj^pi^ 1 )/ 2 v q = 1 and the entries of u satisfy the equations above. 
Clearly, b j~ pqi - q ~ V) l' 1 = (— l)?(9-i) so the eigenvalues of DC P equal i/jfe = au k ~ x with k = 1,2, and where a = 1 
if p(q — 1) is even and a = e m l q if p(q — 1) is odd. Let A = diag(^i, i>2, Vq) and let E be the g x q matrix whose 
columns are the corresponding normalized eigenvectors, computed using Equations I54[ so that DC P E = EA. Then 
Equation 1531 implies that 



p(x)(U ordkr (K, a, A, 0)) = exp [ - i2n G{l,x)] E exp [ - i2 K \ ( & (0 + a/2) (x) A) ] E 1-1 



(55) 



-z2k c(a:) 



-i2«Ac(/3) 



i2rec(x+i) 



£7 _ 



-i2«c(x+S_i) 



2«A C (/3+^_i)) 



where /? = x + 6* + o:/2 + and /i = e l27r< ^. Equations and Inequalities completely analogous to those in Theorem [6] 
can be derived for the on-resonance double kicked rotor operators. 



To illustrate these results, Figures [THE] show examples of the various operator's spectrums. Figure Q] illus- 
trates (j(Uh (k, a, A)) , Figure [2] illustrates a{UkH (k, a, A)) , and Figure [3] illustrates <r(U or dkr(K, a, A)) , for the 
following parameter values: n = 0.25, 0.5, 1, 2, 4, 8, a = 8/13, and A = 1. A grid of 100 values of both x 
and were used for each k to compute estimates for each of these spectrums. Thus, each estimate consists of a 
subset (having 100 x 100 x 13 = 130, 000 points) of the spectrum. Since the operators are unitary, their spectrums 
are subsets of the unit circle in the complex plane, and we have plotted the estimated spectrums on circles to 
reflect this geometry. However, we plot the estimated spectrums on circles whose radii increase with k to enable 
them to be compared visually. We recall from Section UJ that the spectrum of the almost Mathieu mother operator 
<t(H(k, a, A)) consists of q disjoint intervals if q is odd and consists of q — 1 disjoint intervals if q is even. Since 
Uh(k, a, A) = exp ( — inH(a, A)) , the Spectral Mapping Lemma 1181 implies that <j(Uh{ki a, A)) equals the image 
of the union of these intervals under the map t — > exp(-iKt). Therefore, for sufficiently small values of the 
spectrum (t(Uh(k, a, A)) also consists of q disjoint intervals if q is odd and consists of q — 1 disjoint intervals if q 
is even and the lengths of each of these intervals is proportional to |«|. As \k\ increases, some of these intervals 
will begin to overlap and eventually each will wrap around the unit circle and the spectrum will equal the entire 
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unit circle. These tendencies are clearly illustrated in Figure [T] A comparison of Figures [2] and [3] shows that 
the spectrums of the mother kicked operators (for identical parameter values) are identical, thereby illustrating 
Theorem [5j Moreover, it can be shown, using the Baker-Campbcll-Hausdorff formula and operator splitting, that the 
Hausdorff distance between these spectrums and <j(Uh (k, a, A)) decreases with 0(k 3 ), which can be assessed from 
comparing the figures. Figure [4] illustrates <7([/ff(ft, a, A, 9)) , Figure [5] illustrates a{UkH(^, ol, A, 0)) , and Figured] 
illustrates a(U or d,kr{^i ot, A, 6)), for the following parameter values: K — 0.25, 0.5, 1, 2, 4, 8, a = 8/13, A = 1 
and fixed 9 = 0. In these Figures 0HH1 10, 000 values of x were used to calculate each estimate for the spectrums, 
yielding the same total number of points in each spectrum estimate as in the other figures. For a fixed value of 9 1 
the spectrums are proper subsets of the spectrums of their respective mother operators. Even though the spectrums 
for their associated mother operators are identical, the spectrums of o^C/fc/f (k, a, A, #)) and o~(U or< ikr{^, ct, A, #)) for 
a fixed value of 9 are different, as is apparent by comparing Figures [5] and [6] 



IV. FUTURE RESEARCH 

One class of future research projects is to improve existing algorithms and to develop new algorithms for computing 
spectrums. For rational a, a direct method that estimates the spectral bands by directly computing their endpoints 
would be more efficient than our current method that forms unions of sets of eigenvalues of matrices over a grid of 
values of x and 9. Moreover, new algorithms for computing global topological properties of the eigenvalues, such as 
the Chern numbers computed by Leboeuf et al. in [34|, are desirable. 

Another class of projects will address the physical significance of spectral characteristics. These include con- 
ductivity, diffusion and other transport phenomena. 

Finally, a longer term project will address the conjecture that for irrational a the spectrums of all of these 
operators are Cantor sets (regardless of the value of n). Preliminary numerical results supporting this conjecture 
are presented in Figures [THHJ Figure [7] shows the combined width, W, of all the spectral bands of UkH{n,a, A) 
as a function of q with fixed k = 1 and for several values of A. In particular, it seems that W(q) for A = 1 
approximately follows a power law in q with negative exponent, suggesting that the spectrum's measure could be 
zero for A = 1 and irrational a. Figure [5] shows the imaginary part of the logarithm of the spectrum of Ukn{^, ct, A) 
for k = A = 1. This way of displaying the data was chosen for easy comparison when zooming in on different parts 
of the spectrum, and simply amounts to plotting on a line rather than the circles used in Figures HH6] The figure 
uses a — p/q — 2584/4181, which approximates the reciprocal (v5 — l)/2 of the Golden Mean (>/E + l)/2, and can 
be obtained from truncating its continued fraction expansion. The columns show how one obtains strikingly similar 
structures by repeatedly zooming in on the central region of the spectrum. 

We briefly outline a project strategy for proving the conjecture of the spectrums being Cantor sets. The proof that 
the spectrums of almost Mathieu operators for irrational a are Cantor sets is based on the fact that these operators are 
represented on £ 2 (Z), via the Fourier transform, by tridiagonal (Jacobi) matrices. Therefore, the Fourier transforms 
of the generalized eigenfunctions that correspond to the points in the spectrums are sequences in £°°(Z) that satisfy a 
three term recursion relationship that is expressed using 2x2 transfer matrices. Unfortunately, the kicked operators 
can not be represented, via the Fourier transform, by tridiagonal matrices or even by banded matrices. However, a 
powerful result about loop groups can be used to approximate the kicked operators by unitary operators that can be 
represented on £ 2 (Z) by banded matrices and hence enable the use of transfer matrices. This result is Proposition 3.5.3 



in the book Loop Groups by Pressley and Segal [43j which states "If G is semisimple, then L po iG is dense in LG. 



Here G means a semisimple group of matrices, such as the group SU(q) of q x q unitary matrices with determinant 
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one, LG is the loop group of G that consists of continuous functions from T into G under pointwise multiplication, 
and L po iG denotes the subgroup of LG consisting of matrix- valued functions each of whose entries is a trigonometric 
polynomial. This result can be understood as a version of the Weierstrass approximation theorem that asserts that 
every function in C(T) can be uniformly approximated by trigonometric polynomials. This result was used by Lawton 
in [33( to derive a result about a class of filters that can be used to construct orthonormal bases of wavelets having 
compact support, and used in combination with operator splitting methods by Oswald in 35] to show that the error 
of approximating of an element a £ LG by an element b G L po iG decreases with it smoothness of a and the degree of 
b in roughly the same manner as for the Weierstrass's theorem. In particular, since the loop group elements related 
to kicked operators are real analytic functions from T into SU {q) , the approximation errors decay exponentially with 
the degree of the approximating loop group elements in L po [ SU(q). This degree corresponds to the size of the transfer 
matrices. We note that loop groups have already found extensive applications in physics ranging from the Toda 
equation \^ to M-theory 
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Appendix A: PHYSICAL CONSIDERATIONS AND EXPERIMENTAL REALIZATIONS 

This appendix reviews some of the physical considerations and difficulties in experimentally realizing the operators 
considered in Sections [Til an d IIIII 

The first family of operators treated above, H(a, A, 9) in Equationfl] arises from the Harper model, which is a tight- 
binding model describing electrons on a two dimensional square lattice with lattice spacing a, subjected to a uniform 



magnetic field with magnitude B and direction perpendicular to the lattice; see 29] for a detailed description. The two- 
dimensional problem can be separated, yielding Equation [TJ in one direction, and plane- wave behavior perpendicular 
hereto. The parameter 9 is proportional to the momentum of this plane-wave, and the parameter a is given by 

a = ~J~~B, (Al) 

where e is the elementary charge and h is Planck's constant. For parameters describing realistic materials, a « 2A, so 
to have a ~ 1 would require B w 10 5 T - orders of magnitude larger than what can be obtained in present laboratories. 
Nevertheless, other physical situations provide realizations of the Harper model with non-vanishing a, yielding the 
first experimental indication of a Hofstadter's butterfly spectrum in [451 ]. 

Following the rules of quantum mechanics, the Unitary Harper operators Un(n,a, X,9) defined in Equation [3] can 
be seen as the time propagator for a system with a time-independent Hamiltonian given by H(a, A, 9), where the time 
interval propagated is proportional to n. Thus, since the Harper operators have been experimentally realized, so have 
these unitary operators. 

In contrast, the unitary kicked Harper operators [/&#(«;, a, A, 9) defined in Equation [4] have not yet been experimen- 
tally realized, even though they play an important role in the field of quantum chaos. The possibility of experimentally 
realizing operators with the same spectrum as the kicked Harper operators is what currently makes the on-resonance 
double kicked rotor particularly interesting from a physical point of view. 
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The unitary on-resonance double kicked operators defined in Equation [5j 



Uordkr( K i a : \ &) = ex P — «2k cos(27rx) T(a) exp — i2k A cos (2n(x + 9)) T(a) , 

were proposed to be experimentally realized in a cold- atom setup 1511 , Apart from the choice of experimental pa- 
rameters, similar experiments have been performed earlier, see e.g. [30(. Specifically, the proposal considered a long 
quasi-one dimensional gas, consisting of cold atoms or a condensate, subjected to a time-periodic sequence of potential 
kicks caused by a standing wave generated by two counter-propagating pulsed laser beams. The kicking potential is 
sinusoidal with spatial period equal to half the lasers wavelength, and x signifies position within such a spatial period. 
The time dependent Hamiltonian is thus composed of brief time intervals where the dynamics are dominated by the 
kicking potential, and in between these kicks, time intervals where the dynamics are given by free (kinetic energy 
only) propagations. This Hamiltonian commutes at all times with translation by half a wavelength, so assuming that 
the initial state also has this periodicity, the state at all later times will have it as well. We can thus restrict our 
attention to one period x £ [0, 1). We shall have more to say of this periodicity assumption below. 

The operator U or dkr(K, ct, A, 9) is the propagator for one time-period of kicks and free propagation, say from time 
tj to time tj + \. Reading the terms in U or dkr{K, a, A, 9) from right to left, the j'th time-period starts at time t = tj 
and contains a free propagation for a time interval r cx a followed by a first kick at time t = tj + r. Hereafter, a 
second free propagation follows for a time interval Td — t, and a second kick at time Td ends the propagation at time 
tj + \ = tj + Td- The delay time Td will be chosen to be the so-called resonance- or recurrence-time explained below. 
We turn now to look at the individual terms in U or dkr{^, ot, A, 9) in more detail. 

The kicks are the cause of the terms exp [— i2n cos(27ra;)] and exp [ — i2n A cos [2-k{x + 9))], where 9 is the spatial 
separation between the maxima of the two kicking potentials. This separation can be experimentally realized by phase 
shifting the laser pulses. 

The two terms T(a) and T(a)^ 1 are due to free propagation between the kicks. Whereas T(a)~ 1 appears di- 
rectly so in the propagator, T(a) actually arises from the free propagation for a time interval Td — t, yielding 

Therefore, choosing the delay time equal to the resonance time, Td = the first term 

which can be ignored since it equals 1 when applied to any basis function £ n (x). 



exp 



iT d d 2 
2irT R Ix 1 



exp 



— ia d 
47T dx 2 



2ir dx 



in the operator becomes exp 

In passing, we notice that the existence of such a time is due to the free propagation having eigen-energies, which are 
all an integer multiple of a certain energy. We also note that the experimentally realizable ranges for the parameters 



have been estimated to be a £ [0.005, 2], n £ [0.1, 100] and kA £ [0.1, 100], whereas 9 can take on all values [52 1. 

Returning to the assumption that the initial state is periodic in x, we now consider a somewhat more general 
initial state, namely a Bloch wave ^(x) = ij){x) exp(«27r 9b x), where the wavefunction *5f(x) on the entire real axis 
is a product of a function ip(x) having period 1, and a complex exponential with quasimomentum proportional to 
9b- This quasimomentum in the wavefunction would correspond to a real momentum of the entire gas or condensate 
cloud in the experimental setup. As above, the translational symmetry of the Hamiltonian ensures that 9b does not 
change in time. Again, we can restrict our considerations to functions on the interval [0,1) with the basis {£ n (x)}, 
provided we make the substitution d/ dx — » dj dx + i2n 9b- Making these substitutions, and again choosing Td = Tr, 
we see that 



U or dkr(K, a, A, 9) 



exp 



i2it9l 



U, 



ordkr 



\K) a, A, + 6b{cx, — 2)^ exp 







) exp 


'-2 9 B f 

dx 



Thus, the introduction of 9b gives rise to three changes. First, an extra numerical term is introduced, which has the 
effect of moving the eigenvalues on the circle. Second, in U or dkr(n, a, A, 6), the parameter 9 — > 9 + 9s{a — 2). Third, 
an extra operator exp [— 2 9b d/dx] is right -multiplied on U or dkr- With this last extra term, the resulting operator 
no longer belongs to the rotation C*-algebra, and its spectrum can be very different from the original operator 55 1. 
Notably, although the 9 in the kicked Harper operator, Equation [H can be thought of as coming from introducing 
a quasimomentum by d/dx — * d/dx + z27r(9/a, this procedure has little connection to the quasimomentum in the 
experiment proposed in 51 1. 
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An interesting question for future research is thus how this assumption of Ob — influences the experimental 
realizability of U or dkr(n, a, A, 0), and if it can be relaxed. In a real experiment, the cold atoms or condensate is 
trapped in a weak external potential, making the initial wavefunction a wide gaussian rather than a Bloch wave, 
further complicating the situation. That being said, in the usual experimental situation, the wavefunction's finite 
width gives negligible contributions to the outcome of many presently performed measurements. 

Finally, it should be noted that if one were instead to use a ring-shaped trap, the quasimomentum is exactly zero, 
and realizing the operator U or dkr{i^, ct, A, 9) would not be troubled by the issues raised above. Rather than using kicks 
from a laser, one must then use potential kicks that are linear in space, which would yield exactly the cosines in the 
angular coordinate of the ring, and the 0-parameter would be the angle between the directions of the spatially linear 
kicks. 



Appendix B: SPECTRAL THEORY AND C*— ALGEBRAS 



This appendix summarizes results in 3l| and 37| that explain the relationship between the structure of an algebra 



and the spectrums of elements in it. It also derives auxiliary results used in Section |TT] and Appendix [Cl 

Definition 7 An algebra is a complex vector space A equipped with a bilinear map A 2 A, (A,B) -> AB, (called 
multiplication) that satisfies A(BC) = (AB)C. An algebra A is called abelian if AB = BA, unital if it has an (unit 
or identity) element with Ij^A = AI_& = A and A G A is called invertible if there exists an element B such that 
AB = BA = Then B is unique, is called the inverse of A, and is denoted by A^ 1 . The set of all invertible elements 
is denoted by Inv(.4), and for every element A G A its spectrum is defined by 

a A (A) = {^G C : nI A -A & Inv(A)}, (Bl) 

and its spectral radius is defined by r(A) = max{ : fj, S o~(A) }. A subspace J C A is called a subalgebra if 
A, B € J =>• AB G J. We define the center of an algebra A by 

Z(A) = { A G A : AB = BA for every B eA}. 

Clearly 2(A) is a subalgebra of A. A subset J C A is an ideal of an algebra A if A G A, B G J =>■ AB,BA € J. 
The ideals {0} and A are called trivial and an algebra with only trivial ideals is called simple. A homomorphism from 
an algebra A to an algebra B is a linear map <f> : A — > B that satisfies <j)(AB) — (p(A)<j)(B). <f> is called unital if A and 
B are unital and 4>(Ia) = Its- The kernel ker(</>) = { A G A : 4>(A) = 0} is an ideal and <f> is injective if and only if 
kei((f>) — {0}. <f> is called an isomorphism if it is bijective, and then we say that A and B are isomorphic and write 
A = B. If J C A is an ideal then the set of co sets A j J is an algebra under the multiplication (A + J)(B + J) = AB + J 
and the map p : A — * A / J defined by p(A) = A + J is a homomorphism and kcr(p) = J. Every homomorphism 
4> : A — > B induces an isomorphism (p : A /ker(</>) — * <p(A). 

Lemma 8 If A is a unital algebra and B is a subalgebra that contains 7^, then o~a(B) C a^(B) for every B G B. 
Proof. If pb $ <jb(B) then pbl^ — B = pI B - B G Inv(^4) hence p £ a^(B). ■ 

Definition 9 An involution on an algebra A is a map A — > A* that satisfies (A + B)* = A* + B* , (AB)* = 
B* A* , A** — A, (cA)* = cA* and I* A = I a whenever A is unital. Then A is called a *-algebra. If A and B are 
*-algebras then a homomorphism cf> : A —>■ B is called a * -homomorphism if it satisfies <p(A*) = (0(A))*. An element 
A of a ^-algebra A is called normal if AA* = A* A, self-adjoint if A* = A, and unitary if A G Inv(yt) and A* = A^ 1 . 

Definition 10 A normed algebra is an algebra together with a vector norm \\ ■ \\ that satisfies ||AB|| < \\A\\ \\B\\. In 
this paper we only consider continuous homomorphisms between normed algebras. A homomorphism that preserves 
norms is called an isometry. A Banach algebra is a complete normed algebra. If A is a unital Banach algebra and 
A G A then C(A) denotes the (necessarily abelian) norm completion of the set of all polynomials in and A and is 
called the Banach algebra generated by A. 
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Definition 11 A Banach *-algebra is a Banach algebra with an involution * that satisfies \\A*\\ — \\A\\. A C* -algebra 
is a Banach ^-algebra that satisfies = ||^4|| 2 - If A is a unital C* -algebra and S C A we define the C* -algebra 

generated by S to be the smallest C* -subalgebra of A that contains I a o,nd S. If A is a unital C* -algebra and A 6 A 
we letC*(A) denote the C* -algebra generated by {A}. Clearly, C*(A) is the (necessarily abelian) norm completion of 
the set of polynomials in IaiA and A*. If A and 23 are unital C* -algebras, then Hom(_4, 23) will denote the Banach 
space, i.e. the complete normed vector space, of * -homomorphisms from A to B equipped with the operator norm 
topology. If A and B are unital then we require that homomorphisms map I a to I& . 

Examples. C(T), L°°(T), and l°°(Z) are Banach spaces under the norm || • Hoc, i 1 (T) and £ l (Z) are Banach spaces 
under the norm || • ||i, and L 2 (T), £ 2 (Z) and C 9 , q e N are Hilbert spaces. For n £ Z we define S n : Z — > { 0, 1 } by 
5 n (m) = 1 for m — n and S n (m) = for m ^ n. The Banach spaces i 1 (T) and ^(Z) are Banach *-algebras whose 
multiplication is convolution and whose involution is defined by f*(x) — f(—x). Only ^ ,1 (Z) is unital and neither are 
C*-algebras. For instance the function / = -26 -\ + 8 Q + 5i G t(Z) satisfies ||/ * /*|| x = 12 < 16 = ||/||?. The 
Banach spaces C(T), L°°(T), and £°°(Z) are abelian C*~algebras whose multiplication is pointwise multiplication 
and whose involution is complex conjugation. The Hardy subspace H°°(T) C L°°(T), consisting of functions / such 
that (£„ | /) = whenever n < 0, is a Banach subalgebra of L°°(T), but not a *-subalgebra. The algebra 23(H) of 
bounded operators on a Hilbert space H whose multiplication is composition, whose involution is the adjoint map 
and whose norm is the operator norm, is a unital C*-algebra. It is abelian if and only if dim(H) = 1. It is simple 
if and only if q = dim(H) < oo and then it is isomorphic to the algebra A4 q of q by q matrices. The set /C(H) 



of compact operators is an ideal and a simple C*-algebra ( 37[, Example 3.2.2) and every automorphism of /C(H) 
is an inner automorphism, viz. of the form A — > UAU^ 1 for some U € £> M (H), see ((sB], Theorem 2.4.8) and ([3], 
Lemma V.6.1). Furthermore, /C(H) is unital if and only if q = dim(H) < oo and then /C(H) is isomorphic to M. q . 
The C*-algebra C(H) = 23(H) / /C(H) is called the Calkin algebra and it is simple. If dim(H) = oo then this algebra 
has unitary elements that have no logarithm ([37j|. Example 1.4.4). 

Lemma 12 If A is a unital Banach algebra and A € A, then o~a{A) is a nonempty subset o/D(||A||V 



Proof. Lemma 1.2.4 and Theorem 1.2.5 in 37 1. 



Lemma 13 (Beurling) If A is a unital Banach algebra, then the spectral radius of every element A £ A satisfies 

r(A) = inf ;p»||* = lim \\A"\\±. (B2) 

n>l n— >oo 



Proof. Theorem 1.2.7 in [37]. 



Lemma 14 If A is a unital Banach algebra and 23 C A is a Banach subalgebra that contains I a, then for every 
B £ 23, cta(B) C ctb(B) and das{B) C daA(B) (where the boundary dX = closure(X) D closure(C\X) for X C C). 
Furthermore, if o~a(B) has no holes (viz. its complement is topologically connected) then cja(B) — <T]s(B). 



Proof. Lemma [5] implies that o~a(B) C a^(B). Theorem 1.2.8 in 37[ implies the remaining assertions 



Example. If B = z e B = H°°(T) cA= L°°(T), then a A (B) = T c C D(l) = <r B (B) and 9D(1) = T c = dT c . 

Definition 15 (Gelfand) The spectrum (as defined on page 15 in \31 1 ) of an abelian Banach algebra A is the set 
fl(A) of nonzero continuous homomorphisms from A to C equipped with the weak* topology (explained in the Appendix 
in fSj] ). The Gelfand transform of A S A (also defined on page 15 in fSTl l ) is the function A : Q(A) — > C defined 



by A(uj) = lu(A), w G O. 

Lemma 16 (Gelfand) If A is an abelian unital Banach algebra, then 
1. Q(A) is nonempty and compact, 
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2. A £ C(Cl(A)) and a (A) = A(fl(A)) for every A £ A, 

3. the map A — > A is a norm-decreasing homomorphism from A into C(Vl(A)). 

Furthermore, if A € A and if A is generated by I A and A, then A is abelian and the map A : £l(A) —> o~(A) is a 
homeomorphism {i.e. a continuous bijective map whose inverse is also continuous). 

Proof. Theorems 1.3.3, 1.3.5, 1.3.6 and 1.3.7 in 0. ■ 

Example. The Gclfand transform for ^ 1 (Z) is the Fourier transform T : ^(Z) — > C(T), which is defined 
by T(f) — X^riez f( n )£n- It is injective. However, since the function / = — 2<5_i + Sq + 6\ £ ^ 1 (Z) satisfies 
ll-? r (/)l|oo = \/35/2 < 4 = ||/||i, it is not an isometry. J r (£ 1 (Z)) C C(T) consists of continuous functions whose 



Fourier series are absolutely convergent. Wiener 54| first proved that if h S JF(^ 1 (Z)) and h never vanishes then 

l/he 

Lemma 17 (Gelfand) If A is a non-zero abelian unital C* -algebra, then the map A — > A is an isometric 
^-isomorphism from A onto C(D,(A)j, If A is a unital C* -algebra (not necessarily abelian) then r(A) — \\A\\ for 
every normal element A. If A is a unital C* -algebra (not necessarily abelian) and B C A is a C* -subalgebra that 
contains I a, then o~ A (B) = o~b(B) for every B £ B. 

Proof. Theorem 2.1.10 in 1371 implies the first assertion directly and implies the second assertion by applying it to 
C*(A). Theorem 2.1.11 in [37j is the third assertion. ■ 

In view of Lemma 1171 we will denote the spectrum of an element A of a unital C*-algebra simply by o~(A) 
since the spectrum is independent of the algebra with respect to which it is defined by Equation IB1I Furthermore, 
we observe that if A is a unital C* -algebra and if A S A is normal, then every £ tt(C(A)) can be extended to 
4> € Cl(C*(A)) by setting <f>(A*) = 4>(A). The map <j) — ► (j) is a bijection from ft(C(A)) onto ft(C*(A)). Therefore, we 
may identify il(C(A)^j with Q(C*(A)^ and then Lemma 1 161 implies that A may be identified with a homeomorphism 
of fl(C*(A)) onto a (A). Then we may identify A with the restriction of the identity map z : C — > C to a(A) and 
identify A* with the restriction of ~z to cr(A). Therefore, we may identify C*(A) with the C*-subalgebra C*(z) of 
C(a(A)) that is generated by the functions l,z and z. The Stone- Weierstrass Theorem implies that the algebra of 
functions on a (A) that are defined by polynomials in z and z is dense in C(a(A)) and hence that C*(z) — C[a(A)j. 
Therefore, we may identify C*(A) with C(a(A)). 

Lemma 18 (Spectral Mapping Lemma) If A is a unital C* -algebra, A E A is normal, and f £ C(o~(A)), then 
f(A) e A and a(f(A)) = f(a(A)). 

Proof. As noted above, since A is normal C*(A) is an abelian C*-algebra that can be identified with C(a(A)) and 
A can be identified with the restriction of the identity function z to cr(A). Therefore, f(A) can be identified with the 



element f(z) £ C(a(A)) C A. The second assertion is Theorem 2.1.14 in 37j. ■ 

Proposition 19 If A is a unital C* -algebra then d[a(A), cr(B)\ < \\ A — B || for normal A,P> E A. 

Proof. Let s = d[a(A), <t(B)\ and assume to the contrary that s > \\ A— B ||. Without loss of generality (by swapping 
A and B) we can assume that there exists an a £ o~(A) such that s = min{ \a — b\ : b £ cr(B) }. We observe that the 
spectrum of B is disjoint from the open disc of radius s centered at a. Therefore, the operator a/4 — B is invertible and 
normal and its spectrum is disjoint from the open disc of radius s centered at 0. Furthermore, we observe that the map 
z — > z^ 1 maps the union of {00} with the complement of this open disc onto the closed disc D(l/s). Therefore, we 
can apply Lemma IT51 to obtain a((alj\ — B)^ 1 ) C D(l/s). Lemma [TTl implies that IKa/4 — -B) _1 || = r((al^ — B)^ 1 ) 
and hence ||(a//i — -B) _1 || < 1/s. Using the multiplicative property of the norm given in Definition [TO] we obtain 

\\(aI A - A)(aI A - B)- 1 - I A \\ = \\ (aI A - B)' 1 ((aI A - A) - (aI A - B)) || < 1 1 (aI A - B)- 1 1 1 \\A - B\\ < 1. 
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A standard result ([37(, Theorem 1.2.2) implies that (oIa — A)(aIj L — B) 1 £ Jnu(.A) and hence a/4 — A £ Jnw(^4). 
This contradicts the fact that a £ ct(j4) and concludes the proof. ■ 

Proposition 20 Let A and B be united C* -algebras, let X be a compact topological space, and let p : X — > Hom(.4, B) 
be a continuous map such that for every A £ A there exists an x £ X such that p(x)(A) 7^ (then we say that the 
subset p(X) C Hom(yl, B) separates points). Then for every normal A £ A, 

a(A)= |J a(p(x)(A)). (B3) 

xex 

Proof. Let A £ A be normal. If /x G C is not in the spectrum of A then pi a — A has an inverse and hence 
for every x £ X, pig — p(x)(A) has the inverse p(x) [(pi a ~ an d therefore p a(p(x)(A)). This proves that 

a(p(x)(A)) C cr(A) for every x £ X. We will now show that the opposite inclusion holds. We first observe that since A 
is normal, Lemma fTTl implies that without loss of generality we may assume that A — C*(A) . For every idwe denote 
p(x)(A) by A x and the subalgebra p(x)(A) C B by A x . Clearly A x = C*(A X ). We observe that each homomorphism 
p(x) : A — > Ax induces a continuous map p(x) : il(A x ) — > fi(^l) defined by p(x)(<f)) — <j> o p(x), <j> £ f2(„4 x )- 
The discussion following Lemma IT7l shows that we may identify VL(A) with a(A), A with C(cr(A)), and each A x 
with C(<7(A,;)). With these identifications p(x) : cr(A x ) — > cr(^4) is the inclusion map. Furthermore, Lemma IT51 
implies that for each / £ C^A)) and each x £ X we may identify p(x)(f) with the restriction of / to the subset 
a(A x ) of cr(A). Let F denote the closure of (J xeX a(A x ). Clearly Y C a(A). We now show that p(X) separates the 
points in C(a(A)) if and only if Y = cr(A). If Y 7^ then there exists a nonzero function / £ C^A)) whose 

restriction f\y = 0. Then for every x £ X the function p(x)(f) = f\ a (A x ) = and hence p(X) does not separate 
the points in C(a(A)). To prove the converse assume that Y — a(A) and that / is a nonzero function in C(a(A)). 
Therefore, {J xeX a (A x ) is a dense subset of cr(A) and hence the restriction of / to [J xeX o-(A x ) is not equal to 0. 
Therefore, there exists an x £ X such that p(x)(f) ^ and hence p(X) separates points in C(er(A)). We conclude 
the proof by proving that Y = [\ xeX a(A x ). It suffices to show that [J xeX o-(A x ) is closed. Let p £ Y and let 
Pi £ [J xe x a (Ax), i £ N be a sequence that converges to p. It suffices to show that p £ \J xeX cr(A x ). For every i £ N 
there exists an Xi £ X such that pi £ o~(A Xi ). Since X is compact, we may assume without loss of generality (by 
considering a subsequence of Xi) that there exists a y £ X such that Xi converges to y. Proposition [TO1 implies that 
the map x — * cr(A x ) is continuous. Therefore, the sequence a(A Xi ) of points in the hyperspace Ti(a(A)^ converges 
to a(A y ). Since pi £ a(A Xi ) and d[a(A Xi ),a(A y )\ — > 0, it follows that there exists a sequence yi £ cr(A y ) such that 
\Pi — Vi\ — * 0. Therefore, \yi — p\ — > and hence p = y £ cr(A y ). This concludes the proof. ■ 

Definition 21 A representation of a C* -algebra A on a Hilbert space H is a * -homomorphism p £ Hom(.4, B(H)J . 
A representation p is faithful if the kernel of p equals (the zero subspace o/H). A representation p is irreducible if 
and H are the only closed subspaces of H that are invariant under p(A) for every A £ A. 

Lemma 22 Every C* -algebra is isomorphic to a C* -subalgebra of B(H) for some Hilbert space H, or equivalently, 
if it admits a faithful representation. 

Proof. This is Theorem 3.4.1 in [37|. Such a faithful representation is called a Gelfand-Naimark-Segal representation. 



Appendix C: ROTATION C™ -ALGEBRAS AND THEIR REPRESENTATIONS 



This appendix summarizes results in [9| and [18| that explain the structure of rotation C* -algebras. It also derives 
auxiliary results used in Section HT1 
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Definition 23 A rotation C* -algebra is a unital C* -algebra A that is generated by unitary elements U and V that 
satisfy the commutation relation 



UV = exp(i2Tra)VU 

for some a £ [0, 1). The pair (U, V) is said to be a frame with parameter a for A. 
Example 1. Let q £ N and define u> = cxp(i2ir / q) and matrices F,C,D £ M. q by 



(CI) 



F = 



1 



1 



1 



1 



1 UJ l 



1 u 



1 LO- 2 LJ~ 3 



C = 



1 
1 










1 



1 




D = 



1 

w 
lo 2 







uji- 1 



(C2) 



Then (C, D) is a frame with parameter 1/q for the rotation C*-algebra M. q . Furthermore, if a = p/q £ (0, 1) with p 
and q coprime, then (C P ,D) and (C,D P ) are frames with parameter a for .Mg. The matrix F is called the discrete 
Fourier transform matrix. We observe that CF = FD hence C = FDF^ 1 is diagonalized by F. This property is 
extremely useful for developing efficient algorithms. 

Example 2. Let M : L°°(T) -> B(L 2 (T)) and i? : T -> S(i 2 (T)), be defined by Equations US] and [HI let 
a G [0, 1), and let B a denote the C*-subalgebra of B(L 2 (T)) generated by R(a) and M(£ x ). Then (iJ(a), M(^i)) is 
a frame with parameter a for the rotation C*-algebra B a . 

Proposition 24 For a G [0, 1) there exists a unique {up to isomorphism) universal rotation C* -algebra A a that has 
the following properties: 

1. there exists a frame (U a , V a ) with parameter a for A a , 

2. if (U,V) is a frame with parameter a for a rotation C* -algebra C*(U, V) then there exists a unique 
*-homomorphism tt G Hom(„4 Q , C*(U, V)) such that n(U a ) — U and Tr(V a ) = V. 

Furthermore, this homomorphism is a surjection. 



Proof. A proof based on the Gelfand-Naimark-Segal construction is given in Chapter VI of [181 ] . 



For 5 £ N define the equivalence relation = q on T 2 by (z' 1 ,z' 2 ) = q (^1,22) if and only if there exists (wx,W2) G T 2 
such that w\ = w\ = 1, z[ = w% Zi, and z' 2 — W1Z2. 

Lemma 25 For a = p/q where q G N and p and q are coprime integers and (zi, Z2) G T 2 , let ir Zl .z 2 £ Hom(^4 a , M. q ) 
denote the representation of A a on C q that satisfies 7T Zl)Z2 (J7 a ) = z\ C p and n ZuZ2 (V a ) — Z2D. We observe that the 
existence and uniqueness of it ZltZ2 is ensured by Proposition \ 24\ Then every 7r Zl . Z2 is an irreducible representation, 
every irreducible representation is unitarily equivalent to some 7r ZliZ2 , and 7iy <z > and tt Zi . Z2 are unitary equivalent if 
and only if (z[, z' 2 ) = 9 (zi,z 2 ). 

Proof. This is Theorem 1.9 in ■ 

Lemma 26 Let A be a rotation C* -algebra with a frame (U,V) with parameter a = p/q. Let tt G Hom(yl Q ,yl) be 
the homomorphism that satisfies ir(U a ) — U and n(V a ) — V. Then the center 2(A) = C*(U q , V q ) and A is a finitely 
generated Z(A)-module with basis { U r V s : 0<r, s < q — l}. This means that for every A £ A, there exist a unique 
set of elements A r , s £ Z(A), < r, s < q — 1 such that 

9-1 

A= A r , s U r V S . 
r.s=0 

Moreover, tt is an isomorphism if and only if Z(A) is isomorphic to C(T 2 ). This is equivalent, to the condition that 
n(Z(A)) = T 2 . 
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Proof. The first two assertions are proved in the discussion following Theorem 1.10 in [9]. The third assertion is 
Proposition 1.11 in 9]. The last assertion follows from the identifications explained in the discussion that follows 
Lemma [T7] ■ 

Let A be a rotation C*-algebra with a frame (U,V) with parameter a — p/q. For (21,22) £ Tj?, let 
7T z 1 ,z 2 € Hom(A a ,Mq) denote the homomorphisms defined in Lemma 1251 and let tt £ Hom(A a ,A) denote the 
homomorphism defined in Lemma [21)1 Let tt denote the restriction of tt to Z(A a ) and let tt Zi ,z 2 denote the restriction 
of % ZUZ2 to Z(A a ). Lemma [TT] implies that Z{A a ) = C(Q,(Z(A a ))) and Z{A) = C(Q{Z{A))). Since Lemma 
implies that Q(Z(A a )) — T 2 it follows that Z(A a ) = C(T^). We can use the surjection tt to identify Q(Z(A)j with 
its image in T 2 under the injective map <f> — > <fi o tt, £ £l(Z(A)), and to identify tt with the map that restricts 
functions in C(T 2 ) to the closed subset f2(Z(_4)). Likewise, Z(M q ) — C, r2(Z(7W 9 )) consists of a single point, and 
Z(M q )=C(Sl(Z(M q ))). Wc may use the surjection tt Zi ,z 2 to identify fi(.Z(.M,j)) with the subset { (21, 22) } C T 2 . 

Proposition 27 Let ,4, 7r 7 and tt Zi , Z2 be as above and make also all the identifications as above. Then there exists a 
homomorphism p Zl ,z 2 G Hom(.4., M. q ) that makes the following diagram commute, i.e. tt Zi ,z 2 = Pz lt z 2 t, 

M q 

if and only if (21,22) € ^(2(-4)). i/w's case p Zl>Z2 is o,n irreducible representation of A. Furthermore, every irre- 
ducible representation of A is unitarily equivalent to some p Zl .z 2 and the representations p z ' ,z' 2 an d Pz t .z 2 are unitary 
equivalent if and only if (z[, z' 2 ) — q (21, 22). 

Proof. For (21, 22) £ T 2 there exists a p Zl ,z 2 such that ir ZltZ2 = p Zl .z 2 ° 7r if and only if 

kernel(7r) C kernel(7r zi) z 2 ). (C3) 
Lemma [26] shows that Inclusion IC3I holds if and only if the following inclusion holds 

kernel(7r) C kernel(7r ZljZ2 ). (C4) 

Since the map tt restricts functions in C(T 2 ) to the subset Sl(Z(A)) and the map tt Zi ,z 2 restricts functions in C(T 2 ) 
to the subset { (21, z 2 ) } it follows that Inclusion IC4I holds if and only if (z 1? 22) £ Q,(Z(A)). The remaining assertions 
follow from Lemmas 1251 I 

We consider the Hilbert space L 2 (T 2 ) with the standard orthonormal basis 

{€m,n(x>y) = £m(x)€n(y) : {x,y) £ T 2 , m,n £ Z} . 

Then the group SL(2, Z) (integer matrices with determinant one) admits the representation $ : SL(2, Z) — > 
B(L 2 (T 2 )) given by 



rn.n — ^am-\-bn,crn-\-dnj 



a b 
c d 



GSX(2,Z), m,neZ. (C5) 



This representation is called the Brenken-Watatani automorphic representation of SL(2, Z) on A a - 
Lemma 28 A a admits a faithful representation on _L 2 (T 2 ) such that 

U a £m,n = exp(?7ran) £ m +i,n , (C6) 
V a £ m ,„ = exp(-i7ram)^ m ,„ + i. (C7) 
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Furthermore, if <f>(g) and g £ <SX(2,Z) are described by Equation I C5\ then 

U a ^(g)- 1 = exp(~nraac)U^V^, (C8) 
$(.g) _1 = exp(z7ra6d)^t/ Q d . (C9) 

Proof. The commutation relation follows from the computation 

£4 V a Cm,n = exp(-i7ram) f/ Q £ m , n +i = cxp(-i7ram) exp [i7ra(n + 1)] £ m+ i,„+i = exp(i27ra) V a U a £ m>n . (CIO) 

The unitary equivalences follow from almost identical calculations. The first one is 

$(.g) U a $(g)~ l £m ,n = exp [ma(-cm + an)] $(g) €dm-bn+i,-cm+cm = exp(-iiraac) V££ m ,n- (CH) 

■ 

Proposition 29 Define W a = exp(— iira) U a V a . There exists an L £ B(L 2 (T 2 )) such that 

LUaL- 1 = W a , (C12) 
LV a L~ l = V a . (C13) 

Proof. Follows from Lemma [2"51 by choosing L — &(g) where g has entries & = 0, a = c = d = 1. I 

Remark Choosing a = d = 0, b = —1, c = 1 gives the Fourier automorphism &(g)(U a ) — V a and $(g)(V a ) = U^ 1 
that was used in Corollary 2.2 in Q to derive the Aubry- Andre 0] identity a(H{a, A)) = A a((H(a, A -1 )) where 
H(a,X) is defined in Equation [26] 

Lemma 30 For every 9 £ [0, 1), there exists a unique C* -algebra homomorphism ng : A a — * B a such that ng(U a ) = 
ti(ff)R(a) andir 8 (V a )=M(£ 1 ). 

Proof. Since (£i(6)R(a), M(£i)J is a frame with parameter a for B a , the assertion follows from Proposition[24l ■ 

Lemma 31 For 6 £ [0, 1) the homomorphism ng is injective if and only if a is irrational. If a is irrational, then irg 
is a bijection, A a = B a , and therefore 

a(A) = a(Trg(A)), A £ A a , 8 £ [0, 1). (C14) 

Proof. Theorem 1.10 in 0] implies that if a is irrational then A a is simple and hence ng is injective. If a is rational 
then wg is not injective since iigiU^ — £i(qQ) I a J) — 0- The second assertion follows immediately. ■ 

Proposition 32 If fi > 0, < ax < a 2 < 1, and (Uj, Vj) is a frame in A aj for j — 1, 2, then there exists a Hilbert 
space H and injective homomorphisms TTj : A aj — * B(H) for j = 1,2 such that 

27 

\Ui(Ux) - k 2 {U 2 )\\ < 9A*V2^(a 2 -ai), IMVq) - tts(V^)|| < — y/2it{a 2 - ax). (C15) 

A* 



26] 



Proof. Corollary 3.5 in [9J] gives a proof based on the perturbation result ([9j, Theorem 3.1) which is due to Haagerup 
and R0rdam 



Definition 33 A C* -algebra A is called an Azumaya C* -algebra if it is isomorphic to the algebra of continuous 



sections of a matrix bundle whose base space equals il(Z(A)), see This means that there exists an integer q £ N 

and a vector bundle with total space E(A) and projection map 4> : E(A) — > fl(Z(A)) such that for every w £ Q,[Z(A)\ 
the fiber <j> (w) C E(A) is isomorphic to the algebra Ai q . Here, each element A £ A is identified with a function 
A : fl(Z(A)) — > E(A) such that cj> o A is the identity map on £l(Z(A)). With this identification, if A, B £ A and 
w £ n(Z(A)), thenA*(w) = (A(w)), (AB)(w) = A(w)B(w), and (A + B)(w) = A(w) + B{w). 
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Proposition 34 Assume that geN and that A is an Azumaya C* -algebra whose fibers are isomorphic to M. q . Then 
the irreducible representations of A have the form A — > A{w) for w E Q{Z(A)) and where A is identified with a 
section of the bundle described in Definition ] 3 S\ If A G A then 



a(A) = |J eig(AM). 



(C16) 



If fl(Z(A)) is connected then for every A E A, o~{A) is a union of at most q disjoint connected subsets of C, and if 
A is unitary then o~(A) consists of at most q disjoint arcs (or bands) in the circle T c . 



Proof. The bundle representation for rational rotation C*-algebras was constructed in [28|, elaborated in [11], and 
used to obtain a structure theorem for rational rotation algebras in [li^ . Let I q E A4 q denote the identity matrix. 
The second assertion follows since if A E A is identified with a section of the bundle, then for every A E C, A — XI A 
is invertible if and only if A(w) — XI q is invertible for every w E n(Z{A)). The second assertion also follows from 
Proposition [20] The last assertion follows since for every A E A, the set eig [A(w)] is a continuous function of 
wEVL(Z(A)). ■ 



Azumaya algebras were introduced by Azumaya Q (see also 20 1) and C*-Azumaya algebras are discussed in 
0], 41 1- Concrete C* -algebras of operators were used (informally) by Werner Heisenberg and (more formally) 
by Pascual Jordan in the 1920's to model the algebra of observables in quantum mechanics. A special class of 
C*-operator algebras, viz. those closed in the weak operator norm topology, were studied by John von Neumann 
and Francis Murray in a series of papers between 1929-1949. Abstract C*-algebras were invented in 1947 by Irving 
Segal [46J based on the ideas in the seminal 1947 paper by Israel Gelfand and Mark Naimark 23|. A C*-algebra 
A is called almost finite (AF) if it contains an increasing sequence A n , n E N of C* -subalgebras whose union is 



dense. These algebras were introduced in the seminal paper by Bratteli [10]. Rotation C*-algebras appear implicitly 
in Peierls' 1933 paper ([39], Equations 48-53). Irrational rotation C* -algebras were systematically studied by Rieffel 



44j . In 1980 Pimsner and Voiculescu 42j showed that for irrational a, A a can be embedded in the AF algebra 



constructed from the continued fraction expansion of a and in 1993 Elliot and Evans 2l| determined the precise 
structure of irrational rotation algebras. Applications of rotation algebras (also called noncommutative tori) to the 
Kroneckcr Foliation and to the Quantum Hall Effect are discussed by Connes 17]. 
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Figure 1: Spectrums for mother unitary Harper operators, Uh(k, a, A), for a — 8/13 and A = 1. As stated in detail at the end 
of Section [IIII each spectrum is a subset of the unit circle. The figure plots several of these spectrums on concentric circles 
whose radii increase with k G {0.25,0.5,1,2,4,8}. The real and imaginary axes are indicated at the common center of the 
circles. These estimated spectrums are computed using a grid of 100 points in both x and 6, using Equation 1491 




Figure 2: Spectrums for mother kicked Harper operators, UkH (k, a, A), using the same parameters as in Figure [T] and using 
Equation [50] for computations, also with the same (x, #)-grid as in Figure[T] Comparing the spectrums in this figure with those 
in Figure [3] clearly illustrates Theorem [5] namely that the spectrums are identical. 
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Figure 3: Spectrums for the mother on-resonance double kicked rotor operators, U or dkr(K, a, A), using the same parameters 
as in Figure [T] and the same (x, #)-grid for computations. Theorem [S] says that the spectrums in this figure are identical with 
those in Figured] 




Figure 4: Spectrums for unitary Harper operators, Uh (k, a, A, 9), for a — 8/13, A = 1,0 — and for several values of 
k £ {0.25, 0.5, 1, 2, 4, 8}. As stated in detail at the end of Section IIIII each spectrum is a subset of the unit circle. The figure 
plots several of these spectrums on concentric circles whose radii increase with «. The real and imaginary axes are indicated at 
the common center of the circles. We note that these spectrums are proper subsets of the physically more relevant spectrums 
displayed in Figure [T] These estimated spectrums are computed using a grid of 10, 000 points in x, using Equation 1451 
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Figure 5: Spectrums for kicked Harper operators, UkH(n, ot, X,8), using the same parameters as in Figure [4] and using Equa- 
tion 2U] for computations, also with the same i-grid as in Figure [3] The spectrums shown here are proper subsets of those 
shown in Figure [2] Compare these spectrums with the spectrums for the on-resonance double kicked Harper operators, also 
for 9 = 0, shown in Figure HJ] 




Figure 6: Spectrums for the on-resonance double kicked rotor operators, U or dkr(K> Oi, A, 6), using the same parameters as in 
Figure [4] The spectrums shown here are proper subsets of those shown in Figure [3] In contrast to the spectrums of the 
mother on-resonance double kicked rotor operators shown in Figure [3] which are identical to those of the mother kicked Harper 
operators shown in Figure [2] these spectrums differ from the spectrums of the kicked Harper operators for fixed 8 = shown 
in Figure [5] 
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Figure 7: The figure shows a log-log plot of the combined width, W , of all spectral bands as a function of q for three different 
values of A = 2/3, 1 and 1.2. In these calculations, k = 1 and the error for each data point is smaller than 10~ 9 . Interestingly, 
the combined width for A = 1 seems to approximately follow a power-law in q. A fit to these data-points yields the fitting 
parameters 5.26±0.66 for the pre-exponential factor and — 1.22±0.02 for the exponent. This result suggests that the spectrum's 
measure for A = 1 and irrational a could be zero. 




Figure 8: The figures shows e, the imaginary part of the logarithm of the spectrum of UkH(n,a, A) for k = A = 1 and 
a = p/q — 2584/4181. This value of a is a ratio of Fibonacci numbers approximating the reciprocal — l)/2 of the Golden 
Mean (y/E+ 1) /2, obtained by truncating its continued fraction expansion. The first column shows the whole spectrum, and the 
vertical axis indicates its values of e. The second column is a zoom-in of the first column's center part, with the zoom-in factor 
Si = 290. The third to the fifth columns are the zoom-ins of the column preceding them with zoom-in factors of 5*2 = 13.3 
and S3 — Si = 14. There seems to be considerable self-similarity, especially between the second to fourth column. The bands 
of the fifth column are wider than their predecessors, due to the finiteness of q, which indicates that we have reached the limit 
of finding self-similarity by zooming in. 



